HIP-2007-73/TH 



The Partition Function 

bo ; of a Multi-Component Coulomb Gas on a Circle 

O 
(N 

^ '. Niko Jokela,fl Matti Jarvinen,il|and Esko Keski-Vakkuri^ fl 



^Helsinki Institute of Physics and ^Department of Physical Sciences 
P.O.Box 64, FIN-OOOI4 University of Helsinki, Finland 
^ University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark 

Abstract 



O 

♦"^ , We study a two-dimensional Coulomb gas consisting of a mixture of particles 

\ carrying various positive multiple integer charges, confined on a unit circle. 

\ We consider the system in the canonical and grand canonical ensembles, and 
attempt to calculate the partition functions analytically, using Toeplitz and 

^ ' confluent Vandermonde determinants. Just like in the simple one-component 

I ' system (Dyson gas), the partition functions simplify at special temperature 

^ ■ /? = 2, allowing us to find compact expressions for them 

o ■ 

(N 
> 

m 
o 

O 

> 

•i-H 

X 

J3 



* niko . j okela@helsinki . fi 
tmjarvine@ifk.sdu.dk 

* esko . keski-vakkuri@helsinki.fi 



1 



1 Introduction and Summary 



Classical two-dimensional Coulomb gas at finite temperature is a prototype example 
for learning about phase transitions. There are various interesting simplified versions 
of the system. One can confine the particles to move on a unit circle and yet obtain 
an important toy model for various analytical techniques, which also turns out to be 
related to other rather different physical systems pQ. Typically one restricts all charges 
to have the same strength (which can be normalized to unity, leaving the temperature 
as the only parameter) but possibly different signs. Such systems have been solved 
exactly in the bulk and on the circle in the whole stability range of temperatures 
[21 El in isj^ . In the simplest system (the Dyson gas [7j) all the charges are confined 
on the circle and have the same sign and strength. Even the Dyson gas at finite 
temperature is very interesting; among other things it gives a physical realization 
of different random matrix ensembles, just by adjusting the temperature [8]. The 
canonical partition function can be calculated analytically at generic temperature 
and has a simple form [71 [H] , but the grand canonical partition function is much more 
complicated [H [3l |2]. However, at special inverse temperature /3 = 2 it simplifies 
drastically into a very simple form. 

In this paper we consider a slightly more general version of the system, consist- 
ing of a mixture of particles carrying different integer charges, albeit with the same 
(positive) sigii^. This generalization makes the calculations much more complicated. 
At special inverse temperature {3 = 2 there are again simplifications, but now even 
the canonical partition function is more difficult to calculate. We develop techniques 
using confluent Vandermonde determinants, Toeplitz determinants and generating 
functional methods, to arrive at a compact expression for the partition function (at 
(3 = 2) for any combination of po^iti™ integ;^ charge! For the grand canonical 
partition function, we also find a compact expression, but it is more complicated as 
it still contains an infinite series for which we have not been able to find an explicit 
summation. 

Our analysis is mostly restricted to the special inverse temperature (3 = 2 and 
to positive integer charges. It would be important to find generalizations of our 
techniques for all temperatures. This would also allow a treatment of the case where 
some of the particles carry negative charges [H El [3l [10] . Another restriction of the 
present analysis is the number of particles. Although our formulas are exact, they can 
in practice only be used for small numbers of multiply charged particles due to the 
rapidly increasing complexity of the expressions. It is perhaps possible to overcome 
this restriction; instead of exact formulae one may seek asymptotic approximations 

^Note also that a two-component system of charge ratio 1:2 was considered in [6j. 
^Equivalently, one could consider a mixture of Dyson gases, at different (initial) temperatures. 
•^An explicit analytic form can then be found for any given charge expression with the help of a 
straightforward computer algorithm. 
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for a large number of particles, e.g. with the help of asymptotics of large Toeplitz 
determinant^ Together with a generalization to arbitrary temperatures, this might 
shed light on the thermodynamic limit of these systems, including phase transitions 
and analytic properties of the partition functions. 

The article is organized as follows. We introduce the multi-component Coulomb 
gas in the next section. In section 3, we construct the system from the one-component 
Dyson gas by placing several unit charges at the same point and removing the in- 
finite self-energy, and present a result for the canonical partition function of the 
multi-component gas. Technical details of the calculation involving manipulations of 
confluent Vandermonde determinants are postponed to appendix |Al In section 4 we 
develop the former construction into a generating functional method to compute the 
canonical and grand canonical partition functions. We end by a giving an explicit 
example of applying our methods in section 5. 



2 Setup 

The Hamiltonian of the standard Dyson gas with particles is given by 

N 

i/D= 5^ V^(t.,t,) , (1) 

l<i<i 

where 

y(t„t,) = -log|e**'-e**^| (2) 

is the two-dimensional Coulombic potential between the particles i and j which are 
located at e**' and at e**-' on the unit circle, respectively. The definition can be 
naturally extended for particles having multiple positive integer charges. Let us take 
Ni particles of charge +1, N2 particles of charge +2, . . . , N^^^^ particles of charge 
+n^a., labeled by positions {t^^\ . . . , t^J} . • • • , t^}, • • • , • • • , fcj}. 

respectively. They interact via 

Vnm{tf\tf^) = -nmlogle^*'"' -e^^'-'l . (3) 

The Hamiltonian of this system is then given by 

H = Y.^nm{tt\tf^) (4) 



pairs 



N„ 



- ^ ^ log|e**^" - e**-* |- ^ nm ^^lo. 

n=l l<*<i l<n<m i=l j=l 



e ' — e 3 



*We will discuss the behavior of the partition function in the limit of a large number of +1 
charges, in a different framework (decaying branes in string theory) in a future publication 
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The canonical partition function is defined by 



Zc(iVi,iV2,...,iV„_ 



T-max Nn ^^(") 



nn 

n=l i=l 



27r 



(5) 



In this paper we analyze Zq at a fixed inverse temperature J3 = 2. The results are 
also used to calculate the grand canonical partition functio 



•e J 



E i^n^ay^'-- Zc{N,, . . . , N^^J , (6) 



Ni=0 



where Zi are the fugacities which correspond to the different charges. We will start 
by a rather direct evaluation of the integrals over tj in ([5]), and then go on to study 
a more advanced formulation involving generating functionals. 



3 Zq and confluent Vandermonde determinants 

Let us study the partition function ^ at the particular inverse temperature (3 = 2. 
We denote Zi = exp(itj). The partition function of the standard Dyson gas with the 
Hamiltonian ([1]) becomes 



i=l l<i<j 



„ N 



A^! / J-J- 2niz 



'Alzi,...,.^)!^ . (7) 



Here the integrand pd = exp [—2H-^ is the absolute value squared of the Vander- 
monde determinant 



|A(2;i, ...,z^) 



n 

\<i<i<N 



Zi Zj I 



N 



/ J ^«i---2jv^l 
2 



E(-i)"n 



where 11 denotes permutations of 1, 2, . . . , and (—1)^ is the sign of the permutation 



n. Using |Ap = A A* we can express the integrand as an analytic function of Zi 

N 



\A{zi, ...,zn) 



niW-naW 



(9) 



ni,n2 



i=l 



and doing the integrals in ([7]) simply picks the constant term in such that Z-£){N) 
1. 



^For an attempt to calculate the grand canonical partition function of this system, but in a 
different framework, see [12) . 
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The above analysis generalizes to the multi-component Coulomb gas. In fact, 

(10) 

pairs 

can be expressed as a determinant of a confluent Vandermonde matrix (see, e.g., 
|13]). We present here some details of a simple proof since in the case of the Dyson 
gas it has an intuitive physical interpretation. The general treatment can be found 
in appendix A. 




Figure 1: (The square root of) the Boltzmann weight of a state with one doubly- 
charged particle is obtained by forcing two unit charges to coincide. The infinite 
self-energy of the formed pair is removed by a derivative wrt. the angle e between 
the particles. 



The idea of the derivation goes as follows. We start from a single-component 
Dyson gas with at a lot of particles. We then form multiply-charged particles by 
placing several unit charges at the same point and removing the infinite self-energy 
of the formed particle. Let us start from the simplest case of only one double charge 
to see how this works (see also figured]). We take a system with + 2 unit charges 
and denote = tAr+j, Wi = z^j^i. Setting r2 = ri + e we immediately get 

|A(Z1, . . . ,Z7V+2)| = JJ \Zi 

l<i<j<N 

= n = 

l<j<j<Ar 

= [p,{z„...,z^,w,)Y^\ + 0{e') , (11) 

where the leading result gives exactly the wanted Boltzmann weight p2 times the 
exponential of the self-energy that we need to divide away. Notably, as depicted in 



N 

- Zj\ Y\ \Zi - Wi\\Zi - W2\\Wi - W2I 

i=l 
N 

- zj\Y[\zi - wi\h + O {e^) 



i=l 
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figure d] the above procedure equals differentiation, 



-\A{zi,...,ZN+2)\ + Oie) 
e 

A(zi, . . .,Zn,Wi,W2 



W2 - Wi 



+ 0{e) 



d 



dw' 



Zn,Wi,W2] 



(12) 



W2=Wl 



N 



n i=i 

where ([8]) was used in the last step. In particular, from f|T2l) one sees that p can be 
expressed as a determinant of a certain confluent Vandermonde matrix. 



/ 1 



A 



Zl Z2 



1 1 

Zn Wi 



\4 



N+l N+1 
^2 





1 

2wi 



\ 



w 



N+l 



(13) 



(iV + l)< / 



which is obtained from the standard Vandermonde matrix by differentiation. A 
lengthy calculation gives the formula 



ZciN, 1) 



A^!l! 
1 

m 



N 



i=l 



dti dr 



n^^^rP2(ti, 



N 



n 



dzi dwi 



N 



i=i * ^ ni,n2 



-i)"^(-i)"^n 



HiW-naW 



X [ni(Ar + 2) - 1] [n2(Ar + 2) - 1]^ 



ni(7V+l)+ni{Ar+2)-n2(Ar+l)-n2(Ar+2) 



— (iV + 3)(Ar + 2)2(Ar + l) 



4 



+ 



iV + 3 
4 



(14) 



where the complex integrals again pick the constant term of the integrand. The result 
( IT4l) equals the one obtained by using the Selberg integral [12] . 

We present the straightforward, but technical generalization to any number of 
charges in the appendix El In particular, particles carrying higher charges in > 
2) are obtained by taking higher derivatives. The final result for Zq in a general 
configuration is given in fH6|) . Notice that the result can be generalized to some 
special higher values of (3 by adjusting the charge configuration of the particles. This 
is possible if \J l3/2qi are positive integers for all the charges qi of the system. For 
example, Zq for a system with /3 = 8 and for particles of the charge n (with 
n = 1, . . . , rimax) is obtained by setting 



; 



2n 



(15) 
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in (06]). 

Similar results for Zq as those of the appendix |X] can also be derived by using 
a more elegant method that involves a generating functional, which we will discuss 
next. 



4 Generating functional method 



Let us define a generating functional (see appendix [B] for an alternative definition of 

M 

n 

\<i<3 



oo 

M=0 ^ 



i=l 



I 77- of ■ \ 



(16) 



The weight of the standard Dyson gas pd is then found by taking functional derivatives 

of/b1, 

Pd(^i, • • • , ^n) 



(17) 



i=o 



For the special inverse temperature /3 = 2 we can use to write (see, e.g., [H]) 



oo 

= Em/ 



M=0 

oo 



4 = 1 



M 



^(_l)ni(_i)n,-Qg.t.[ni 



(i)-n2W] 



ni,n2 

M 



i=l 



M! 

M=o ni,n2 



Jn2(i)-ni(i) , 



i=l 



where Ilj are permutations of the length M and 

jn = 



Ztt 



(19) 



are the Fourier coefficients of the source. That is, I[j] is a sum over the Toeplitz 
determinants of j. 



/[j] = 5^ det tW[j], 

Af=0 

where T^^^) is an M X M matrix with the entries 



(20) 



(21) 



One might expect that the formula f|T7j) can be generalized to particles with higher 
charges by taking higher functional derivatives at a single point, or letting some of 
the tj's in (flTl) coincide. However, as above in section [3] this leads to a problem with 
the infinite self-energy. Luckily, as above in (1121) . the self-energy can be removed by 
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taking suitable derivatives wrt. tj. For example, the weight for a state with unit 
charges and one double charge can be expressed as 



P2{ti,. . . ,tN,T) = pY){ti,...,tN,T,T + e)e + C (e 



1 

2! 



_d_ 

dT2 



d 



Pi){ti, . . . ,tN,ri, T2) 



r=Ti =T2 



2! \dT2J 6j{h)---6j{tN)6j{n)6j{T2) 



(22) 



i=o 

T=Tl=T2 



where = exp (2i?seif) was added to remove the self-energy. As seen from ( l22l) com- 
bining the differentiation idea of the previous section (fT2l) to the functional derivative 
result f|T7j) leads to a "normalized" second order functional derivative at a single point, 



D2 ~ e 



dt 52 



dr 
2^ 



1 



d 



6' 



2! V^rs; 5jin)SjiT2] 



(23) 



-I T=T1=T2 



where the infinite self-energy has been removed. 

We shall now work out a generalization of (!23|) for D„ with n > 2. It can be 
conveniently expressed in terms of the Fourier modes j„ of the source j{t) [see the 
definition (1281) below]. Let us consider the case where one particle (at e*"^) carries an 
arbitrary integer charge n, which can be added by using the result (|38|) in appendix 
Ri For n > 2 fl22|) generalizes to 



p„(ti, . . .,tN,r) 



n 



fc=2 



2fc-2 



(2A;-2)! V^T-, 

_l^n{n-l)/2 " / (5 ^ 



7^ — ) Pd(^1, • • • , ^7V, Tl, • • • , Tny 



(24) 



n 



PT)iti, . . . , ^AT, Ti, . . . , r„) 



In particular, the definition in terms of derivatives can be written in various forms. 
Any combination is fine as long as the sum of the orders of the derivatives is n{n — 1) 
and the normalization factor K of the appendix |A] can be worked out. We shall use 
the latter form of (1241) . Then the partition function reads 



ZcniN) = ZciN,=N,N2 = 0,...,N^., = 0,N^ = l) 

^_l^n{n-l)/2 



(25) 



N\n\ [{n-iy.y 

n-l 



N 

n 



i=l 



dti dr 
27r 27r 



d_ 



_d_ 



n — 1 



s:N+n 



m 



S3{ti)---SjitN)6j{n)---6j{rn) 



j=Q 

T=Tl = ---=Tn 
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The partition function can now be evaluated by using the result (120|) . The chain rule 
implies 
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m=—oo 



m=— oo 



'imt 



d 



dir. 



(26) 



Only the terms that include products of exactly N + n oi j^'s give non-zero contri- 
bution which fixes M = N + n m fl20l) . Hence 



N\n\ [(n- 1)!] 



dr 



27r ydn 



d 



n-l 



_d_ 



n — 1 



(27) 



mi=—M mn=—M 
N M 



k=l 



mkTk 



i=o 

r=T\ = --- = Tn 



1 



d 



M 



(ELi"^fc,o) 

n! [(n-l)!]" 



.A;=l 



-1 5 



'0/ mi=-M mn=-M 

where 5{m,n) denotes the Kronecker delta. Here the operatoi[^ 



dj 



detT(^-^)[j] 



j=0 



M 



2^ ■■■ 2^ r^irrr^-nir 11^"^^-^ 



<9j 



n\ [(n-l)!]' 

mi=-M m„=-M '-^ ' ^ k=l 

generates the particle with the charge n in the gas. The result fl28p is valid for (3 = 2, 
but depends on temperature, in general. The temperature dependence arises from the 
self-energies which can be removed by differentiation similarly as above at all integer- 
valued inverse temperatures. Thus fl28l) can be easily generalized to all /5 = 1, 2, . . .. 
However, the generating functional can be written in terms of Toeplitz determinants 
as in dlH]) only for /? = 2. 

Generalization of (1271) to any number of multiply-charged particles is straightfor- 
ward. If we set M = Yll2T in the definition ([28]), theifl 



Zc(iVl,...,iVn_J 



n 



n=l 



{Dn) 



N ' 



detT(*^)[j] 



(29) 



i=o 



Inserting the definition ( 1211) one can express Zq as a finite sum similarly as in ap- 
pendix |A] (see also the concrete example below). 

Also, a compact formal expression for the grand canonical partition function ([6]) 
at /3 = 2 immediately follows. 



Zcik, ■ ■ ■ , Zn^aJ = exp 



n=l 



A/=0 



(30) 



j=0 



^Note that the Dn are not unique. For example, using the first fine in (124]) would lead to 

h 

"^Notice that £)i = 



-^n — Z^rni,...,m„ = -j\/ ™fc^O (2fe~2)! a} 
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where the summation hmits in the definition (l28l) need to be extended from ±M to 
±00. 



5 An example 

As a concrete example of the use of ( l29l) continue from (1271) and show how a purely 
combinatorial result is derived in the case of only one multiply-charged particle. By 
inserting the explicit form of det T^^'^"* [j] from (fT8l) we find 



N\ 

M M 



j=0 

M / n \ " 

\n-l 



^ ' mi=-M m„ = -M \k=l / k=l 

N n 

X ^ i-'^)^'i-'^)^^Yl^niii),U2ii)YlSn^ij+n)+m,,n2ij+n) ■ (31) 



Hi, Ha i=l j=l 



In this special case of only one multiply-charged particle, the condition Ylk=i ~ ^ 
follows the others in fl3Tl) . and may be dropped. The conditions = 112(2) for 

i = 1, . . . , N define a certain "subdeterminant" , which may be written as 



[(n-1)!]" ^ ^ 
where 11 is a permutation of length n. 
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M n _^ 

E E(-i)"n[^'^-M"' ' ^^^^ 



A A general combinatorical formula for Zq, 

In this appendix we derive a general expression for Zq using confiuent Vandermonde 
determinants and compute the consequent combinatorical result. For the general 
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multi-component Coulomb gas, ([8]) reads 



2nm 



nkJ"^-^r =|detA|2, (33) 



pairs 



where A is a M x M confluent Vandermonde matrix with M = Yln'^T ^^n- A is 
defined by 

A = (^A(i) (zf)) ■ ■ ■ (^41^ ■ ■ ■ j . . . ) , (34) 

where each M x n subblock A^") corresponds to one particle of the gas and has the 
entries 



(j - 1)! \dzj 1 if i < J 

Equation ([33]) with the definition flMl) . for A is a known result, but for use 
elsewhere we will sketch how it follows from the standard Vandermonde determinant 
(IH]). Let f(ti, . . . , tN, Ti, • • • , T„) be a function such that 

n 

/(ti, . . . , tTV, n, . . . , r„) = ^(ti, . . . , tTV, n, . . . , r„) JJ (xi - Tj)'''^ , (36) 

i<i<i 

where kij are positive integers and g stays finite and non-zero as any — ^ tj. Then 
take the mjth derivative wrt. Tj for each i = 1, 2, . . . , n of the both sides of (136|) such 
that mj > and 



j=l l<i<j 

and let r = Ti = ■ ■ ■ = r„ in the end. On the right hand side only such terms survive, 
where all the derivatives act on the product. Thus 

1 / 9 f d 

g{ti,...,tN,T,...,T) = — i—j •■■(^1 /(ti,...,tiv,ri,...,r„) 



5 

r=ri = ---=T„ 



where is the constant integei^l 



(3^ 



l<i<i k=l 



where K is the coefficient of the monomial r{"^ ■ ■ ■ r™" in rii<i<j(''"« "~ Tj)^^^ . To prove 
( !33|) for the case of unit charges and one particle with charge n, we apply the result 



^The result only makes sense if K ^ 0. 
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dMl) to 

f{Zi,...,ZN,Wi, 



N 



N n 



n ~ ^ n 



= JJ (wi - Wj) X . . . ,2;7v,wi, . . . (40) 

l<i<j 

and take rrik = k — 1. After calculating the integer K = Y[k=ii^ ~ ^V- 



N 



N 



l<i<j 



i=l 



,Zn,W, 



d 



w) 



k-1 



A{zi, ...,Zn,Wi,.. .,Wn) 



n ^ 

± det {zi) ■ ■ ■ A(^) {zn) {w)) 



W = Wl = ---=Wn 

(41) 



where we used the determinant form of A in ([8]) and (!35|) in the last line to complete 
the proof. The general result [equation (1331) with (1341) ] follows by a straightforward 
induction. 

In the rest of this appendix, we use the explicit form flM|) . (153]) of A to find a 
combinatorical formula for Zc, the constant term in all z^"''' of fl33l) when expanded 
into an analytic function of the zl"'\ Note that since we only need the determinant 



of A and \z, 



Hi 



1 we may multiply any column by an arbitrary power of z^"'^ . Using 



in addition elementary column operations A^"^ (z) may be replaced by 



(j - 1)! 



(42) 



Then we find 
IdetAP 



M n n 

{i} 



n=l k=l 



■J--|- y'i(n,k)+s) ^^(„)y^{",fc)+ 



E^-M n n 

{«} n=l k 



1 



where the Vandermonde matrices of the permutation variables in the last form are 
obtained after antisymmetrization. The function 



n-l 



£(n, fc) = E + {k- l)n 



(43) 



m=l 
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is only needed for picking up the permutation variable i with the correct index. 
The constant term of | det = det A* det A is 



• "max I ''max ' n -i 

Zc(iVi, . . . , iV„_) \[N^\= Y1 ^n-^M^n-.M U U rj, TT, ^ (44) 

n=l {f^j} n=lfc=l[^K^-l)!---l!] 

n,k)+li ■ ■ ■ ^ 3l(n,k)+n) 
X5 {il{n,k)+l + ■ ■ ■ + ie.(n,k)+ni jl{n,k)+l + " " " + jl(n,k)+n) 

where = 6ij is the Kronecker 5-symbol. For n = 1 the 6 restrictions give simply 

ik = jk- Using these the result evaluates to 

'^max '^max ^ n -j^ 

Zcll^-'-= E ^-..x^.wxlinn^ 7^2 (45) 

S,th{j} n=2 tl[^K^-l)!---l!] 

X A {S{ii'^n,k)+l), ■ ■ ■ , S{ii'(n,k)+n)) A {S {ji> (^n,k)+l) , ■ ■ ■ , S {] e {n,k)+n)) 
f n n 

(n,fc)+s J 



n=2 



,s=l s=\ 



where K = M — Ni, the first sum goes over all increasing functions 5* : {1, . . . , K} — >■ 
{1, . . . , M} [so that i<] ^ S{i) < S{j)], and i'{n, k) = i{n, k) - Ni. 

Due to symmetry, one can add the restrictions ii'(n,k)+i < 'ie'{n,k+i)+i (for all n > 1 
and 1 < A; < A^„), and i(>'(n,k)+s < ie'{n,k)+s+i, ji'{n,k)+s < je(n,k)+s+i (for all > 1, 
/c, and 1 < s < n) and multiply by the ratio of numbers of terms (Iln=2^ ?t,!^"A^„!). 
Then the result becomes 

^max -j^ 

Zc= E ' ^n-^K^n-^K n n TTTTTT^ ^ (46) 

5,{i},{i} n=2 k=l \\"' -L;- ■ ■ ■ J--J 

X A {S{iii(n,k)+l), ■ ■ ■ , S{iii(n,k)+n)) A {S {] {.n,k)+l) , ■ ■ ■ , S {je> (n,k)+n)) 

(n n 
'{n,k)+s ) 
s=l s=l 

where the prime indicates the presence of the above restrictions. 



B Relation to field theory 

In this appendix we present a compact formulation (for all inverse temperatures /?) of 
the generating functional f|T6|) using the language of field theory. We use the operator 
X{t) that lies at 2; = e** on the boundary of the unit disk. Its self contraction 

(X(ti)X(t2)) = log|e**^-e^*H' (47) 
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is essentially the Coulomb potential between charges at e**^ and at e**^. Then the 
generating functional (1161) reads 



m 



exp 



dt 



— j(t):eV/5/2XW 



as can be verified by using the Wick theorem. The canonical partition function of 
the (single-component) Dyson gas reads 



I ^ : exp (vW^X(t)) 



1 

iV! 



dt 5 



N 



m 



j=0 

(49) 



where the last step can be checked using fj48|) . 
In the case of the multi-component gas 



Zc{N, 



N ' 

1 ' rimax ) 



1 



so that one would naively expect that 
Zo(A'i,...,A'„..,„; ^ 



n 

n=l 



dt 
2^ 



: e 



nJl3/2X{t) 



(50) 



n=l 



n 

n=l 



dt S" 



1 Nn 



m 



(51) 



j=0 



However, (1511) fails in general since the higher order functional derivatives at the 
same point are not compatible with the normal ordering: the right hand side includes 
self contractions between fields at the same point. These are the self-energies of the 
multiply-charged particles which were discussed in the text. Equation flSTl) holds 
only after the subtraction of the self-energies, which was done for (3 = 2 explicitly in 
section m leading to ( l29i) . Notice also that the grand canonical partition function can 
be written as (see [T2] ) 



exp 



n=l 



dt 



Pl2X{t) . 



(52) 



which explains the exponential form of the result fl30l) . 
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